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The complex patterns of the hadronic spectrum have puzzled physicists since the early discovery 
of the "particle zoo" in the 1960s. Today, the properties of these myriad particles are understood 
to be the result of quantum chromodynamics (QCD) with some modification by the electroweak 
interactions. Despite the discovery of this fundamental theory, the description of the hadronic spec- 
trum has long been dominated by phenomenological models, due to the difficulties of addressing 
QCD in the strong-coupling regime, where nonperturbative effects are essential. By making nu- 
merical calculations in discretized spacetime, lattice gauge theory enables the ab initio study of 
many low-energy properties of QCD. Significant efforts are underway internationally to use lattice 
QCD to directly compute properties of ground and excited-state baryons. Detailed knowledge of 
the hadronic spectrum will provide insight into the character of these states beyond what can be 
extracted from models. 

In this review, I will focus on the latest progress in lattice calculations of the Pii(1440), the poorly 
known hyperon spectrum and the energies of highly-excited states of the nucleon. Delta and other 
light-fiavor baryons. In the heavy-fiavor sector, I will concentrate on recent lattice-QCD calculations 
of baryon masses, particularly those that make predictions concerning yet-to-be-discovered baryons, 
such as ficc, or triply-heavy baryons. 

00 



o 



00 

o 
o 



PACS numbers: 12.38.Gc,14.20.Gk,14.20. Jn,14.20.Mr,14.20.Lq 



^ i. introduction 

Oh 

^ Hadron spectroscopy plays an important role in understanding the theory of quantum chromodynamics (QCD), 

which has been successful in describing many aspects of the strong interaction. At the most basic level, the SU(3) 
nature of the color interaction predicts the existence of certain kinds of color-neutral hadrons: quarks and antiquarks 
can bind into mesons, and three quarks can form baryons. There have been many experimental observations of these 
^ simple hadrons with all kinds of flavor content and a variety of excitations. QCD also suggests the possibility of 
bound states formed only by gluons, called glueballs, or more complex hadrons, such as exotic mesons, tetraquarks, 
pentaquarks, dibaryons and hadronic molecules; many hadronic spectroscopy experiments are looking hard to find 
evidence of such states. 

In the regime of simple meson and baryon spectroscopy, there are complicated patterns of experimentally observed 
states, and although many of them are understood by theory, there remains much work to be done. Even for light- 
quark hadrons, made of up, down or strange quarks only, there are many excited states, especially in the baryon 
sector, whose physical properties are poorly understood. One famous example is the nature of the Roper resonance, 
• • A^(1440) Pii, which has been the subject of interest since its discovery in the 1960s. It is quite surprising that the rest 
. energy of the first excited state of the nucleon is less than the ground-state energy of nucleon's negative-parity partner, 

the A^(1535) Sn [1], a phenomenon never observed in meson systems. There are several alternative interpretations 
^ of the Roper state, for example, as the hybrid state that couples predominantly to QCD currents with some gluonic 
^ contribution |2 or as a five-quark (meson-baryon) state [3]. 

There are also baryonic excited states predicted by the quark model, based on three constituent quarks but not yet 
discovered in experiments, the so-called "missing resonances", which also provide challenges for the experimentalists 
and also for the phenomenologists. Many experiments at, for example, Jefferson Lab, MIT-Bates, LEGS, MAMI, 
ELSA, and GRAAL, have been collecting data from one and two-meson production to broaden our understanding of 
excited states. These processes are so complicated that many theoretical centers at, for example, Jefferson Laboratory, 
George Washington University, Bonn, and Mainz, have built models to be able to analyze the experimental data. These 
joint efforts will provide us detailed understand of the QCD theory and paint a more complete picture of hadronic 
resonances. 

In the heavy-quark sector, even many of the ground states, such as doubly bottom or triply charmed baryons are 
yet to be discovered. Experimental and theoretical studies of charmed and bottom hadrons have been the focus of 
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vigorous research over the last several years fl'-'7| , such as recent experimental discoveries of new charmed baryons 
at SELEX O [9]. In the summer of 2007, CDF[6] reported the first observation of the heavy baryons and E^, 
and then both DO [7] and CDF [ID] observed the bottom baryon SjJ", breaking a long period of silence following the 
observation of the A;, in 1991. The bottom-baryon spectrum has become somewhat controversial due to recent 
results from DO and CDF. Last summer, DO^ll^ reported a first observation of the doubly strange bottom baryon 
at 6.165(10)(13) GeV. However, a recent CDF work puts the fl^ mass at 6.0544(68)(9) GeV, a difference of 
111(12)(14) MeV. With a discrepancy of 6.2 standard deviations, it appears that the two collaborations cannot both 
be observing the fib- A theoretical understanding of the bottom baryon spectroscopy from first-principles QCD is 
crucial and can help disentangle such discrepancies in experiment. 

The Beijing Spectrometer (BES-III), a detector at the recently upgraded Beijing Electron Positron Collider 
(BEPCII), has great potential for accumulating large numbers of events to help us understand more about charmed 
hadrons. The antiProton ANnihilation at DArmstadt (PANDA) experiment, a GSI future project, and the LHC is 
also expected to provide new results to help experimentally map out the heavy-baryon sector. For these reasons, 
lattice quantum chromodynamics calculations of the spectrum of heavy baryons are now very timely and will play 
a significant role in providing theoretical first-principles input to the experimental program. At the Tevatron, much 
data remains to be fully analyzed; one could possibly expect more discoveries over the next couple years. Furthermore, 
it is anticipated that in the upcoming dedicated bottom-physics experiment LHCb at CERN or a future Super-B fac- 
tory, there will be many more discoveries of bottom-hadron properties. Combining such experiments with improved 
theoretical understanding, our knowledge about these states will in the near future be significantly enhanced. 

Although significant progress has been made in refining phenomenological models since the development of the 
original quark model in the 1960s, determining the properties of hadrons from the first principles of the Standard 
Model requires us to directly solve QCD in the strong-coupling regime. Unfortunately, when the coupling constant 
becomes large, perturbative techniques based on Feynman diagrams do not work well. Instead, we use a discretization 
of space and time in a finite volume to calculate the observables of QCD numerically, a technique known as lattice 
QCD. 

Lattice quantum chromodynamics (LQCD) has been used successfully to compute many experimentally observable 
quantities from first-principles calculation of Euclidean-time hadron correlation functions, even occasionally predicting 
experimental results before they are measured. However, the early successes of LQCD were mostly restricted to 
computation of the physical properties of the lowest-energy states in each quantum number channel. This is due to 
the focus on the large-time behavior of correlation functions, where uncertainties due to excited-state contributions 
are exponentially suppressed. Given that signal-to-noise in correlation functions also falls exponentially at large times, 
success is often dictated by available computational resources. 

With the most state-of-the-art supercomputers available, we are just beginning to be able to simulate at the 
physical pion mass. A recent work by the BMW collaboration [121 calculating multiple lattice spacings, volumes and 
pion masses as light as 180 MeV provided an excellent demonstration of how ground-state hadron masses with fully 
understood and controlled systematics are consistent with experiment. 

In this review, selected channels of light- and heavy-flavor baryons efforts from lattice QCD calculations arc sum- 
marized and compared. In Sec. |lTj a brief introduction to lattice gauge theory is given, along with some examples 
of operators and correlator construction. A few commonly used methods used to extract masses of hadrons are also 
presented. Lattice-QCD calculations of excited baryon spectroscopy involving light and strange quarks are presented 
in Sec. |III| The status of Roper-resonance studies is summarized. We describe progress of the Hadron Spectrum 
Collaboration (HSC)'s effort on highly (both radially and orbitally) excited states, along with a few remarks on recent 
developments and future efforts. We also cover lattice-QCD calculations in hyperon channels with 2 -|- 1 dynamical 
ensembles. In Sec. |IV[ charm and bottom-baryon spectroscopy from various lattice-QCD groups are summarized. The 
heavy flavor is simulated with relativistic heavy quark action, NRQCD and static quarks with various light-fermion 
actions for the light /strange quarks. In Sec.|vj we summarize the current state of the held and discuss some aspects 
of future developments. 

II. THEORETICAL BACKGROUND 
A. Lattice-QCD Setup 

Lattice QCD is a discrete version of the continuum QCD theory. The path integral over fleld strengths at inflnitely 
many Minkowski space-time points from continuum QCD is approximated using only flnitely many points in a discrete 
Euclidean space lattice with periodic boundary conditions. In momentum-space, this automatically provides both an 
ultraviolet cutoff at the inverse lattice spacing and an infrared cutoff at the inverse box size, and it is important to 
make sure that these cutoffs do not affect the physics. Since the real world is effectively continuous and infinitely 
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large, a lattice calculation needs to take limits where the lattice spacing a — >■ and the simulated spacetime volume 
y — > cx) to eliminate the artifacts introduced in a discretized finite box. 

At a fixing lattice spacing, the systematics due to the finite lattice spacing can be controlled by a systematic 
expansion called Symanzik improvement |13' for the lattice fields, actions and operators. The goal is to construct 
an effective continuum action using operators that are invariant under discrete rotations, parity-reversal and charge- 
conjugation transformations, that represent the long-distance limit of the lattice theory and zero the leading finite-a 
errors. However, the breaking of continuous (Euclidean) SO (4) symmetry allows many new degrees of freedom, leading 
to various lattice actions that return to the same continuum action once the symmetry is restored. Thus, there exist 
many gauge and fermion actions for us to choose from. 

For example, in constructing a relativistic heavy-quark action on the lattice, one starts with the Dirac operator and 
the gluon field tensor (and the presence of the heavy quark mass distinguishes between the time and space components 
of each). There are eight operators with dimensions up to five: 



However, not all of the above operators are relevant. That is, the redundant continuum operators may have their 
coefficients arbitrarily set without affecting the resulting physics. The remaining coefficients of the relevant operators 
will then be tuned to give the correct physics in the continuum limit. 

Today, most gauge actions used are 0(a^ )-improved and leave small discretization effects (©(a'^Agj-.j-,)) due to gauge 
choices. On the other hand, most fermion actions are only 0(a)-improved and have systematic errors of 0(a^ Aqq^) 
that become dominant. For this reason, lattice calculations are generally distinguished according to the fermion action 
used. Differences among the actions are benign once all systematics are included, and the choice of fermion action 
is constrained by limits of computational and human power and by the main physics focus. The commonly used 
actions are: domain-wall fermions (DWF) P^17j . overlap fermions|18|. clover (or Sheikholeslami-Wohlert) act ion [TO], 
twisted- Wilson fermions [5U] and staggered fermions [STi 1^ . 

For a calculation done at a typical fixed lattice spacing, 0.12 fm, the discretization effect for an 0(a)-improved 
fermion action is on the order of 0(AQQj-,a^) « 3%. However, one should keep in mind that it varies with fermion 
actions in the actual calculation and the physical quantities (for example, the charmonium hyperfine splitting has 
stronger lattice-spacing dependence than the proton mass). At a fixed lattice spacing, we may expect different amounts 
of discretization effect for different fermion actions; however, after continuum extrapolation, that is, using multiple 
calculations at different lattice spacings to extrapolate to zero lattice spacing, we expects universal results for all 
quantities. 

In order to make predictions using QCD on the lattice, we calculate observables corresponding to vacuum expecta- 
tion values of operators O, taking the form 



Sg is the gauge action and Sp — ipMip is the fermion action with Dirac operator M . The bilinear structure of the 
fermion action allows the integration over the fermion fields to be done explicitly, bringing down a factor of det M. 
This means that the anticommuting fermion fields (impossible to simulate on a computer) are integrated out, leaving 
an integrand that depends only the values of the gauge fields. In the early days of lattice QCD, the computational 
resources were insufficient to compute the fermionic determinant; instead, the determinant was approximated by a 
constant. This is equivalent to removing quark loops from the Feynman diagrams of a perturbative expansion, and 
this technique became known as the "quenched approximation" . Althought quenching retains many of the important 
properties of QCD, such as asymptotic freedom and confinement, it introduces an uncontrollable systematic error, 
and modern calculations keep at least the up, down and strange quarks in the sea. 

The discrete integral we have derived can be evaluated numerically using Monte Carlo methods. The Monte Carlo 
integration uses random points within the domain of gauge configurations to approximately evaluate the integral. 
The "importance sampling" technique is introduced to perform this task more efficiently: instead of choosing points 
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FIG. 1: A selection of time-slices showing the topological charge of the vacuum calculated on one gauge-field configuration 
on a lattice with spacing approximately 0.1227 x 10~^^ meters with 16 lattice sites in each spatial direction. The evolution 
of the vacuum shown in this image occurs over 5.8 X 10~^^ seconds. The red regions correspond to instantons with positive 
topological winding number, while the blue correspond to instantons with negative winding number. 



from a uniform distribution, they are chosen from a distribution proportional to the Boltzmann factor e~^°'^'--'^\ which 
concentrates the points where the function being integrated is large. 

Using this method, we accumulate an ensemble of gauge field configurations, generated using a Markov-chain 
technique. Based on the current state of the gauge configuration, a new configuration is selected. A transition 
probability P{[U'] -s— [U]) is determined based solely on this new configuration and the current one. The new 
configuration is added to the ensemble or rejected, according to the resulting probability. If the "detailed balance" 
condition 

Pi[U'] ^ [C/])e-^^«(^) = P{[U] ^ [[/']) e-^"«(^'' (3) 

is satisfied by this update procedure, then the canonical ensemble is a fixed point of the transition probability matrix. 
Under this condition, repeated updating steps will bring the gauge-field distribution to the canonical ensemble. 

Figure [l] shows a few timeslices of the spatial distribution of the topological charge of a QCD vacuum ensemble 
taking from one of the lattices generated by HSC. 

Once we have accumulated gauge-field configurations using this procedure, the expectation value of our observable 
(O) in the original field theory is numerically approximated by the sample average on our ensemble: 

1 ^ 

(o)^j^J2oiu;:{x)). (4) 

n 

Some of the primary observables we study on the lattice are two-point correlation functions. A particle of interest 
is created by an operator carrying the appropriate quantum numbers. For example, meson states can be created with 
local bilinear operators as shown in Table [ij Baryons may be constructed in the same way; for example, for nucleon 
spin-1/2, we can write down 3 local-site (that is all 3 quarks resides within the same lattice point) operators: 

x'^^Hx) = e-''^{u^^'^{x)Cd\x))j5u'{x) 

X^^\x) = e'''''{u^^''{x)C-f5lid'ix))u-{x) (5) 
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TABLE I: Meson states created by local operators of the form ipFip, labeled in spectroscopic notation 
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FIG. 2: A selection of timeslices showing the real part of the proton two-point correlation function calculated on one gauge-field 
configuration on a lattice with spacing approximately 0.1227 x 10^^^ meters, having 16 lattice sites in each spatial direction. 
The difference between adjacent timeslices corresponds to At = 1.17x 10~^^ seconds, and the evolution of the proton correlation 
function shown in this image occurs over 7.02 x 10~^* seconds. To enhance detail, the colors (red being more positive and blue 
more negative) have been normalized on each timeslice. The proton source on timeslice t = has a Gaussian distribution in 
space. 

More complicated hadrons (such as pentaquarks) and multi-particle states simply increase the number of quarks and 
antiquarks involved. 

For a meson created at a specific space-time y = (y,0) (point source) and annihilated at a; = (^j^) (point sink), 
the corresponding Euclidean two-point function is 

C'^'\x,y) = {0{x)0{y)) = ^J [dC/][#][#]e-^-(^''^^^)(V^(x)r^(x))(V;(y)r^(2;)). (6) 

Figure [2] shows a few timeslices of the spatial distribution of a proton two-point correlator on a QCD vacuum ensemble 
generated by HSC. 

When t is large enough, the lowest-energy (ground) state dominates, and its mass can be obtained from fitting the 
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correlator to an exponential: 

C^'Ht) = (^X^C(2)((f,i),(y,0))^ ^ Are-™-*. (7) 

The so-called "effective mass" can be obtained from In C''^-' {t) / In C'^^-' {t + 1) ; plots of the effective mass as a function 
of time are helpful in crosschecking the fitted values with different fitting methods and adjusting the fitting range 
away from excited-state contamination. However, such a character makes it difficult to study excited states on the 
lattice, since their signals will become exponentially dominated by lower states at large times; later, I will describe a 
few efforts to overcome this difficulty. 

Very often, we impose a smeared source on the quark field using gauge-invariant smearing or gauge-fixed configu- 
rations to improve overlap of the creation and annihilation operators with the desired states: 

y 

with smearing function F. 

Jackknife analysis often adopted to analyze lattice data. It is is a method of regrouping a data sample (similar 
to bootstrap analysis) in order to make a reliable estimate of the error for a quantity computed from that sample. 
Consider an ensemble of N gauge configurations: Ei,E2, ...,En. The k^^ jackknife result for a quantity A is obtained 
by 



- {A)e,^e, (9) 



and the variance of the jackknife estimators is 



(10) 



where A^'^'f is the average of the jackknife blocks. The jackknife analysis is a better estimator when dealing with 
correlated measurements. 

So where does lattice QCD stand in general? 

Ensemble generation is still costly even with the supercomputers available today and the latest algorithmic im- 
provements. The cost of an ensemble increases with smaller lattice spacing, larger volume and lighter pion mass 
composed by the sea quarks. Since light quark masses greatly increase the cost of an ensemble, lattice calculations 
often use unphysically heavy pions. Currently, more groups are generating pion masses around or below 200 MeV, 
and some have started to work at the physical pion mass. When using heavy pions, to get results that connect with 
the real world, we need to extrapolate the pion mass to its physical value. Since chiral perturbation theory is often 
used, this process is generally called chiral extrapolation. One side benefit of this process is that we can help to fix 
the low-energy parameters of chiral effective theory. 

The left-hand side of Fig. [s] shows typical masses (nucleon and fl) extracted from lattice QCD as functions of m^. 
The different-color points indicate calculations done using different lattice spacings; the dashed line indicates the 
physical pion mass. With a lightest pion mass about 180 MeV and no visible lattice-spacing dependence, the chiral 
and continuum extrapolations are well controlled and the results agree nicely with experiment. The right-hand side 
of Fig. [3] shows more hadron masses extracted by various collaborations with different types of discretization due to 
their choices of lattice fermion actions. The excellent agreement among different calculations and with experiment 
for these simple quantities gives us confidence to work on other calculations that are either difficult to measure or 
not-yet measured in experiment. 



B. Extracting Masses 

A variety of analysis techniques have been applied to extracting the hadron spectrum from correlation functions 
in lattice QCD. In this section, we discuss some of these approaches that have been used for extracting excited 
states, including multiple-correlator least-squares fitting, Bayesian fitting [24H26] . black-box methods [27H30] . and the 
variational method[3T| [32]. 
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FIG. 3: (left) The pion-mass dependence of the nucleon and Q, baryon masses at three lattice spacings: a ~ 0.065, 0.085, 
0.125 fm. Figure taken from Ref. [12]. (right) The low- lying spectrum of light hadrons from dynamical lattice-QCD simulations 
from various collaborations around the world with different choices of fermion actions. The solid and dashed lines are the PDG 
values for the masses and widths, respectively. Figure taken from Ref. |23| . 



1. Multiple-Correlator Fits 

The most widely used method is nonlinear least-squares (NLLS) fitting to a model function, such as a sum over 
two or more exponentials. The brute-force way of extracting excited states is to fit using a form that includes them 
explicitly. Usually, one needs multiple correlators as input to make the fit algorithm converge. With multiple smeared 
correlators, one minimizes the quantity 



x'=y2 '^' — (11) 



s.t 



where the correlators are each fit by the sum of A'' exponentials, having independent amplitudes a but the same masses 
(energies) E. The greater the number of distinguishable input correlators, the easier it will be to extract higher states. 

Unfortunately, since the contribution of each state to the two-point correlation function declines exponentially over 
time with exponent equal to the mass, the excited-state contributions are swiftly dominated by the ground state and 
by noise. Indeed, much work where the ground states are of primary interest take pains to only fit the data at large 
times where the excited states have vanished. And since the excited states will have significant contributions only for 
the first few timeslices beyond the creation operator, little data is available to constrain the fits. Thus, this method 
is restricted to very small numbers of fitted masses. 



2. Bayesian Fitting 

One method for extracting additional masses from correlators is to apply more sophisticated fitting techniques. 
Using information we already know, we can introduce what are known as Bayesian priors to constrain our fits. These 
constraints may improve the precision with which difficult-to-determine masses are found by a fitting method. 

The Bayesian technique applies known information to shift the estimated probabilities of correctness for parameter 
sets being considered by a fitting procedure. This prior information may be as simple as the assertion that all 
physical masses must be positive [Ei > 0) or that each new excited state must be heavier than the states preceding 
it {Ei > Ei-i). We could even attempt to extract states from a correlator in an iterative fashion. Assuming we have 
somehow acquired a ground state for our data (say, by a non-Bayesian fit), we can make the prior assumption that 
the ground state is actually equal to that fitted value within the determined uncertainty. 

Using this prior information, we wish to modify the likelihood estimate for all parameter sets considered by our 
Bayesian fitter. The probability of correctness for a particular parameter set is generally represented by the value. 
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which for a non-Bayesian fit is 

2 _ 1 f yz~ f{{a,E}a„n^l,x^)Y 

>^dof-^^l^(^ J^^ ) > (12) 

where we fit N data points y ± 6y to some form / parametrized by n amphtudes a and masses (energies) E. Rather 
than aUowing the fitter to arbitrarily change the amphtudes and masses, we would like it to take into consideration the 
masses we have already determined. Therefore, we simply add a term that will increase (decrease the probability 
of correctness) if the fitter moves away from these values: 

where we have already determined the first m masses to be i? ± SE. 



3. Black-Box Methods 



Operationally, even simple nonlinear fits can be fraught with difficulty, from establishing the range of Euclidean 
times included in the data set to stabilizing the convergence of minimization algorithms by careful choices of initial 
guesses or temporarily freezing selected fit parameters during the minimization process, all of which require manual 
intervention. 

At such times, we often turn to black-box methods for guidance, because they do not require intervention to 
determine initial guesses and fitting ranges: estimates of correlation functions go into the black box and estimates of 
hadron energies come out. The main detraction of black-box methods is that the produced estimates are expected to 
have larger uncertainties, making them suboptimal relative to least-squares methods. Marrying the two approaches 
can effectively combine the best features of both, leading to a highly-automated analysis programs producing optimal 
estimates of energies. 

The problem to solve is the nonlinear system of equations y = V(a;) a 
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for Xm — exp (p)] and a„i — Am (p) exp [— ioE'm (p)] where y„ — C(p, /:„). V(a;) is known as 2M x M 

rectangular Vandermonde matrix. The M = 1 solution is simple to compute and widely known in the lattice QCD 
literature; it is the effective-mass solution. The M = 2 solution was explicitly constructed in Refs. [371 [IHl 133] ■ One 
can analytically solve up to M = 4 systems. 

However, there are ways to get even higher states. For the general approach for M > 4 (Prony- Yule- Walker method 
or just Prony method [34Vl36j ) and solutions with multiple correlation functions, please refer to Ref. [29] for more 
details. 



4. Variational Method 

The variational method |31l I32j is a powerful tool for extracting multiple excited states using a matrix of two-point 
correlators. 

We can often write down multiple independent operators that describe the same quantum numbers. In the simplest 
case, we can merely adjust the smearing function applied to the quarks in a single-site operator. If we include 
operators that involve different sites, there are more operators we can write down to generate even larger bases for 
the operators. We organize these by their behavior under the elements of the cubic group (the group describing the 
rotational symmetry of the lattice) into irreducible representations (irreps). Since the action of QCD respects the 
symmetry of the cubic group, an operator belonging to one irrep never create states having overlap with other irreps. 
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Using these operators, we construct an r x r correlator matrix, Cij{t), where each element of the matrix is a two- 
point correlator composed from different operators Oi and Oj . Then we consider the generalized eigenvalue problem 

C{t)i; ^ X{t,to)C{to)^, (15) 

where the selection of to depends on the range of validity of our approximation of the correlators by the lowest r 
eigenstates. If to is too large, the highest-lying states will have exponentially decreased too far to have good signal- 
to-noise ratio; if to is too small, many states above the r we can determine will contaminate our extraction. Over 
some intermediate range in ipi we should find consistent results. 

If the eigenvector for this system is and a goes from 1 to r, the correlation matrix can be approximated as 

r 

aj = Y.^'*^7^~'''" (16) 

n=l 

with eigenvalues 

A„(t,io)=e-(*-*«)^'' (17) 

(up to exponential corrections with exponents equal to the mass splitting with the next state not included in the 
approximation) by solving 

C{to)-'/''C{t)C{tor'/^^ = X{t, to)V. (18) 

The resulting eigenvalues A„(t,io)i called the principal correlators, are then further analyzed to extract the energy 
levels. En- Since they have been projected onto pure eigenstates of the Hamiltonian, each principal correlator should 
be fit well by a single exponential. The leading contamination due to higher-lying states is another exponential having 
higher energy; a two-state fit may help to remove this contamination. 

III. LIGHT/STRANGE FLAVOR 
A. Roper Resonance 

In both meson and baryon spectroscopy there are many experimentally observed excited states whose physical 
properties, such as spin and parity, are poorly understood. In such cases, theoretical input from LQCD could help to 
solidify their identification. Among the excited nucleon states (N*), the nature of the Roper resonance, A^(1440) Pu, 
has been the subject of interest since its discovery in the 1960s. The reversal of the mass ordering of the positive- 
parity excited state (the Roper) and the negative-parity ground state (A'^(1535) 6*11) is quite unusual, never being 
seen in mesons. This has led phenomenologists to postulate that the Roper might not simply be a radial excitation 
of the nucleon. Rather, it might be a hybrid state that couples predominantly to QCD currents with some gluonic 
contribution PJ or a five-quark (meson-baryon) state Experiments at Jefferson Laboratory, MIT-Bates, LEGS, 
Mainz, Bonn, GRAAL, and Spring-8 offer new opportunities to understand in detail how nucleon resonance properties 
emerge from nonperturbative aspects of QCD. 

Early LQCD calculations using the quenched approximation 1331 137l - t42| , found the computed spectrum inverted 
relative to experiment, with Pu heavier than Sn. Fig. |4] and Table |TT] show a summary of parameters and analysis 
methods used in these works and the extrapolated masses at the physical pion mass. Only Ref. [38 , which uses lighter 
pion masses, seems to observe a potential mass reversal in their central values. There is agreement with experimental 
values within errors, since the masses of the Pu and Sn are overlapping within their statistical errorbar. 

However, if we look at most of the published even-parity LQCD results {Nf = 0) [351 BTHIS] as functions of pion 
mass, we find big discrepancies in the calculated nucleon first-excited mass. The Roper calculation on the lattice is 
in a more chaotic situation than other states. One may suspect that the discrepancy is caused by how the lattice 
spacing is determined. In fact, there are different choices available for the scale-setting technique: Some prefer to 
use experimentally well-measured hadronic quantities, such as f^^ or or bottomonium splittings to set the scale. 
Others prefer to calculate quantities that are not likely to depend on the lattice volume, such as the string tension in 
the heavy-quark potential in the static limit; this is matched up with phenomenological parameters at distance scale 
ro or ri, which requires independent determination. However, if the different ways of setting the scale caused the 
discrepancy we find in Roper masses, we would expect the same amount of discrepancy to be refiected in the nucleon 
mass as well; however, the effects there are much smaller. Refs. [351 US] both use chiral fermion actions (overlap 
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FIG. 4: Summary of previous lattice calculations with extrapolation to the physical pion mass point or the lowest simulated 
pion point (labeled as "f " ) . 



TABLE II: Summary of existing published and Pn calculations. Due to space limitations, we adopt these abbreviations for 
fermion actions; Domain- Wall Fermions [TiHTT] (DWF), Chirally Improved Dirac Operator [JTIHH] (CIDO), Fat-Link Irrelevant 
Clover 49, (FLIC); and for the analysis methods: Variational Method (SU [32] (VM), Constrained Curve Fitting [26| (CCF), 
Maximum Entropy Method [SOI EU (MEM), Black-Box Method gTl [H [52] (BBM). For those works which do not perform 
extrapolation, we use the lightest pion mass to represent their results. 
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fermions and chirally improved fermions respectively), which preserves the chiral symmetry at finite lattice spacing. 
However, their Roper masses are no closer than any other set. Refs. [JSJHS] use local nucleon operators with different 
smearing parameters and analyze the data with the variational method; their results seems to be consistent within 
2(7. Ref. [44] also uses the variational method and adopts more sophisticated operators organized according to their 
representations under the cubic symmetry group. These calculations using Wilson-type fermions are expected to have 
larger discretization error than other 0(a)-improved fermion actions. Taking account the systematics, Refs. [i^Hill^ 
are rather consistent, but disagree with results from Ref. '38^ which adopts Bayesian fitting method. 

We address the same issue on left-hand side of Fig. [5] sorting the results as functions of m^ri, which gives us a 
rough estimate of the expected finite-volume effects on these data. We find that certain discrepancies of the Roper 
masses seems to line up with mj^L. This suggests finite-volume effects can be more severe for excited states than 
ground states and that careful examination of such systematic errors is crucial. However, there are still two regions 
of Roper masses among different lattice groups for m^n-L > 5 and further studies will be required to understand the 
difference between the different ways of extracting Roper masses and properly assign systematics, which are too often 
ignored. 

The examples given above are calculated in a vacuum with giuonic degrees of freedom only; that is, no sea fermion 
loops contribute to the ensembles. In recent years, newer calculations need not apply this approximation, and more 
results with Nf = 2 + 1 (degenerate up and down plus strange in the sea) are available. 

Two preliminary studies were reported in the annual lattice conference using Nj = 2 -|- 1 DWF and clover lattices 
by RBC/UKQCD and xQCD; their results are shown as the open symbols in Fig.l6] but no refined numbers have been 
reported since. CSSM lattice group [54] uses local baryon operators with different smearing parameters to construct 



11 



3.0r 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

) CD 



2.5- 



> 2.0- 
O 



1.0 

0.5 f 


3.0r 

2.5- 



® 



> 2.0- 

(U 

O 
^ 1.5 



1.0 
0.5 








in 4 
4 



I ' I ' I 
4 



B + 

* Mathur (3.2 fm) 

+ Sasaki (3.0 fm) 

O Lasscock (2.6 fm) 

A Burch(2.4fm) 

□ Brommel (2.4 fm) 

O Basak(2.3fm) 

V Mahbub (2.0 fm) 

_] I I I I I I I I 1 I 1 I 1 I 1 I I I I I I I I I I L 



0.2 



0.4 



0.6 



0.8 
M2(GeV2) 



1.0 



1.2 



1.4 



1.6 



a4 ^ 



on 



4 



4 



* Mathur (3.2 fm) 

+ Sasaki (3.0 fm) 

O Lasscock (2.6 fm) 

A Burch (2.4 fm) 

□ Brommel (2.4 fm) 

O Basak (2.3 fm) 

V Mahbub (2.0 fm) 

_l , , , ,_ 

10 



15 



FIG. 5: Summary of published Nf = LQCD calculations of the nucleon and Roper masses in GeV as functions of (left) 
and Roper masses as a function of the dimensionless product of the pion mass and lattice size L (right). Note that the errorbars 
are only statistical. 



a large basis for the nucleon with spin 1/2 and analyze the data using the variational method. They use Nf = 2 + 1 
isotropic lattices at a fixed lattice spacing, a = 0.0907 fm and fixed volume, 32^ x 64 with lightest pion mass 156 MeV. 
Their calculation of the Roper masses are shown in Fig.|6]as 5 purple diamond points. However, the lightest two pion 
masses may suffer finite- volume effects (which have not currently been estimated), since their m^L values are less than 
3. Hadron Spectroscopy Collaboration reported two studies using a large basis of independent operators organized 
according to representations of the cubic group, including operators where the three quarks inside the nucleon are 
located at different lattice sites [55l 156) . They have used 2 + 1-flavor anisotropic clover lattices, with finer temporal 
lattice spacing to improve the excited-state signal in Euclidean time; that is, the lattices have better resolution in 
the direction from which the particle energies are extracted. The green pentagons in Fig. |6] are taken from the lowest 
two states reported in their study last year using only operators in the Gig irrep. This year, a refined study mapped 
the cubic-irrep operators to definite spins using derivative operators constructed with continuum symmetries. They 
found four states closely aligned and roughly around the excited state found by the CSSM lattice group. A refined 
CSSM study using larger basis of their smeared operators reported in this year's Nstar Workshop a similar structure 
as Ref. [SS] . This may suggest that multiple-particle operators need to be included in the correlator matrix to clearly 
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FIG. 6: Summary of the latest Nf = 2 + 1 Roper results as a function of pion mass squared. Note that the filled symbols 
and solid errorbars (open symbols and dashed errorbars) indicate the results taken from published papers (the latest Lattice 
Conference proceedings). 



identify the nature of these closely related states; both HSC and CSSM are currently investigating this. 

The dynamical results in Fig. [6] are impressive but the errorbars shown above are only statistical only; i.e., there 
is no estimation of the systematics due to the finite volume and nonzero lattice spacings in the calculations. This is 
sometimes difficult due to the limited amount of dynamical lattices available. Despite the difficulties, NPLQCD made 
such a calculation in Ref. [57| with high statistics using Nj = 2 + 1 anisotropic clover lattices at ss 390 MeV. Four 
volumes are used in this calculation, 2, 2.5, 3, and 4 fm with TOjri w 3.9, 4.8, 5.8, 7.7. The volume dependence is not 
noticeable in the earlier HSC study with fewer trajectories available at the time. To our surprise, the ground state 
with rriTrL around 4 (rule-of-thumb suggests that at this value finite- volume effects can be ignored) shows around 
30 MeV discrepancy from its value in the infinite- volume limit. They also extract the axial decuplet-octet- meson 
couplings from volume-dependent heavy-baryon chiral perturbation theory, |.gAAf7r| — 2.80(18) (21); similarly for other 
couplings such as |gE*A7r| — 2.49(23)(35), |(7h»hw| = 2.49(23)(35). The effect of the finite volume on the Roper mass 
is even more dramatic; one expects the central values to shift to by about 170 MeV. The multiple-volume study 
also allows the issue of the Roper's multiplicity to be addressed. Since the ratio of the "weight" or overlap factors 
associated with the Roper between the largest two volumes is Afi{4: fm)/Aji{3 fm) = 0.964(52) consistent with one, 
we see that these first-excited states are likely to be one-particle states instead of two. A scattering state would have 
overlap factors proportional to the lattice volume. More time should be devoted to studying volume dependence to 
get a better estimation on the finite-volume effects from these dynamical calculations. 

We can also study the nature of the Roper resonance by examining its form factors, such as the transfer-momentum 
dependence of the Roper-nucleon transition form factors. A first exploratory study was done in 2008, using quenched 
anisotropic lattices at pion mass 720 MeV|58j. Additional pion masses (480 and 1100 MeV) were reported in Ref. |28j . 
However, none the results yielded negative values for F2 *^ in the low-momentum region, as suggested by experimental 
studies. An updated study of the Roper-nucleon form factors using 2 + 1-flavor anisotropic lattices [59l |60j seems 
to reproduce the behavior of CLAS's analysis, providing further strong evidence that the Roper is the first radially 
excited state of nucleon. Further study of the systematics and improved operators will be needed to strengthen this 
statement. 



B. Highly Excited States 



In the previous section, we found how difficult it is to extract the Roper resonance. What if we want to get to 
even higher excited states with LQCD? Unfortunately, in Euclidean space, excited-state contributions to correlation 
functions decay faster than the ground state. Therefore at large times, the signals for excited states are swamped 
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by the signals for lower-energy states. One way resolve this issue is to improve resolution in the temporal direction. 
An anisotropic lattice where the temporal lattice spacing is finer than spatial spacings can provide better resolution 
while avoiding the increase in computational cost associated with a similar reduction of all spacings. 

One also needs large basis of independent operators that allow the maximal overlap with higher-spin operators. The 
use of cubic-group irrep operators or operators constructed according to continuum symmetries is needed to expand 
the reach of the lattice calculation. The variational method becomes quite reliable when allowed to consider up to the 
8th excited state. In this section, we will focus on the Hadron Spectroscopy Collaboration (HSC)'s concerted efforts 
and developments over the past few years. 

1. Dynamical Anisotropic Lattices 

Using a finer temporal lattice spacing helps to reduce the systematics associated with energy extraction for higher 
excited energies (even with the aid of the variational method) . Such lattices provide many time slices of information 
to reconstruct more accurately the signal due to a particular state, which is crucial since the signal in Euclidean time 
declines exponentially. In principle, we could use very fine lattices while keeping the volume big enough to avoid 
the "squeezing" systematic effect; however, computational cost scales significantly (power of 5-6) with the inverse 
lattice spacing. Tuning the parameters for dynamical anisotropic lattices is extremely difficult. Not only do we have 
to set the usual coupling constants and quark masses, but also the anisotropy for both gauge and fermion fields. 
The task is further complicated by the feedback among the parameters. HSC adopted anisotropic clover lattices 
to keep the spatial lattice spacing coarse, avoiding finite-volume systematic error, and to make the temporal lattice 
spacing fine enough to extract towers of nucleon excited states. However, tuning the 0(a)-improved parameters for 
fermion actions in the dynamical gauge generation correctly is a much more difficult task for anisotropic lattices. 
They use Symanzik-improved gauge action and clover fermion action with 3-dimensional stout-link smeared gauge 
fields; the gauge ensemble is generated using the (R)HMC algorithm. The spatial lattice spacing is Os = 0.1227(8) fm 
(determined using too), and the renormalized anisotropy = as/at is 3.5. Ref. |59j gives a detailed study of how 
to set the dynamical anisotropic lattice parameters, and Ref. [60] reports basic lattice properties along with the 
ground-state hadron spectrum. 

HSC tunes the strange-quark mass using the dimensionless parameter sq = ^^^"4^2 "^"^ ■ The parameter sq is set 

as close as possible to the corresponding experimental value at the S'C/(3)/-symmetric point (875-MeV pion mass). 
When the sea-quark mass is reduced, the sji parameter is seen to remain roughly constant for pion masses as low 
as 383 MeV in Ref. [5D]. (Furthermore, the parameter is more sensitive to the strange-quark mass than other 
methods, such as those using the J-parameter. Since it is a ratio, there is no need to worry about estimating the 
shift in the lattice spacing when decreasing the sea-quark mass or generating more statistics.) Ref. |61j reported an 
updated measurement of the same quantity on a 230-MeV 24^^ ensemble, and sq does not deviate from the expected 
value. The sq parameter is a stable and useful observable for setting the strange-quark mass, as shown in the leftmost 
of Fig. |7| 

Ref. ]61j also employs a refined baryon mass extrapolation using too or toh as a reference mass. Using mass ratios 
instead of masses decreases the statistical (and possibly systematic) errors due to explicit cancellation of systematics 
and gauge-field noise and removes the ambiguity associated with setting the lattice spacing. Using the same data, they 
update the mass-ratio chiral extrapolations by modifying the next-to-leading-order heavy-baryon chiral perturbation 
theory. They find the finite-volume corrections are negligible according to ChPT estimates, and the extrapolated 
baryon masses agree with experimental values significantly better than with a linear extrapolation. The results are 
shown in middle and right-hand side of Fig. [7| 

2. Cubic- Group Irrep Operators 

Another important tool for extracting highly excited states is the variational method [3T1 , which uses a matrix 
of different source and sink operators to project more exactly onto the eigenstates of the Hamiltonian and make the 
analysis more reliable. To use the method effectively, we need a large number of independent operators that overlap 
well with excited states with desired quantum numbers and linearly independent of with each other. Unfortunately, 
the single-site baryon operators of form given in Eq. [5]can only give us 3 different operators. We need to expand our 
consideration to a bigger operator space to extract highly excited states reliably. 

The Hadron Spectrum Collaboration (HSC) has been investigating interpolating operators projected into irreducible 
representations (irreps) of the cubic group [44.5 ,62^ in order to better calculate two-point correlators for nucleon spec- 
troscopy. First, we construct various gauge- invariant baryonic operators by displacing (with appropriate gauge parallel 
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FIG. 7: (Left panel) The location of the dynamical ensembles used in Ref. [HT] in the sn-ln plane. The leftmost cir- 
cle (black) indicates the physical point {/n'^'*'^, s^'^'^''}, while the red, green and blue points are the lattices with atms = 
—0.0540, —0.0618, —0.0743 respectively in Ref. [BD]. The horizontal dashed (pink) line indicates constant sn at the physical 
point, and the diagonal line indicates three- flavor degenerate theories. The purple points are from recent measurements of 
the 230 MeV ensemble, showing that it does not deviate away from the Nf — 3 value. Mass-ratio chiral extrapolations as 
functions of the Is and sh for octets (center panel) and In and sn for decuplets (right panel). The lines indicate the "projected" 
leading chiral extrapolation fits in /= and ss {In and sn) while keeping the other one fixed. The black (circular) point is the 
extrapolated point at physical /= and ss (In and sn)- The masses for N, E and A are 0.962(29), 1.203(11), 1.122(16) GeV, 
respectively, and the masses for A, E* and E* are 1.275(38), 1.426(25) and 1.580(13) GeV, respectively. 
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TABLE III: Subduction of continuum baryon spin irreps into representations of the cubic group 



transport) the component quarks in spatial directions as shown in Fig. |8] and compose them with the appropriate 
flavor structure, isospin and strangeness. This large basis of operators does not respect rotational symmetry nor 
the subgroup of rotations that are also respected by the lattice. We apply group-theoretical projections such that 
operators will transform according to an irreducible representation under lattice rotations and reflections. Since the 
action of the lattice theory respects these symmetries, the representation of a state is conserved as it propagates 
through the lattice gauge field; operators belonging to difi'erent representations will not mix even as they move from 
source to sink. Using such a design for our operators, we can create large correlator matrices, and the technique can 
be extended to meson and multi-hadron computations. 

In the cubic group Oh, for baryons, there are four two-dimensional irreps Gig,Giu,G2g, G2U and two four- 
dimensional irreps Hg and I7„. (The subscripts "gr" and "u" indicate positive and negative parity, respectively.) 
Each lattice irrep contains parts of many continuum states. The Gi irrep contains J £ {^,|,|,^,...} states, the H 
irrep contains Js{|,|,|,|,...} states, and the G2 irrep contains J € {1,^,^,...} states. The continuum-limit 
spins J of lattice states must be deduced by examining patterns of degeneracy between the different Oh irreps. For 
example, Gi will be dominated by spin-1/2 baryons with some contribution from highly excited spin-7/2 states. The 
H irrep will contain the spin-3/2 baryons (such as the ground-state Delta), but it will also contain many closely 
spaced states from all higher spins. G2 will contain spin-5/2 and 7/2 states and serves as an important check of these 
spins in the H and Gi channels. See Table [5] for further details about spin and irreps. 

Using these operators, we construct an r x r correlator matrix and extract individual excited-state energies by 
applying the variational method, as described above. More details on the baryon correlation functions evaluated 
using the displaced-quark operators are described in Refs. |44j|62]. Demonstrations of how these operators work using 
purely gluonic vacuum with light valence quarks for nucleon and delta spectroscopy are reported in Refs. 53, 63, 64]. 
Further calculations of isospin-^ excited nucleons in two-flavor QCD, using u and d quarks that have the same mass, 
are reported in Ref. with 2 pion masses: 416(36) and 578(29) MeV. 

One drawback of calculations using these cubic-group operators is that they require multiple orientations in order 
to maximally overlap with a wide range of quantum numbers, and generating so many correlators is quite expensive. 
Furthermore, as we go to lighter and lighter pion masses, there will be increasingly many decay modes open, even for 
the lowest energy at a specific quantum number. We need to extend the matrix to include multiple-particle operators 
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FIG. 8: The orientations of 3 quarks within a single baryon. The right-most one involves all 3 spatial axes. The quark 
displacements shown above use only quarks separated from the central core by one unit; however, larger displacements can be 
used to improve overlap with higher orbitally excited states. 



(so that we can further understand the nature of the "resonance" in our calculation) and "disconnected" operators. 
Further, we need to achieve better precision for each state to distinguish among them. 

A new way to calculate timeslice-to-all propagators, "distillation", has been proposed in Ref. [66 . The method 
is useful for creating complex operators, such as those used in the variational method, since it reduces the amount 
of time needed for contractions and allows the structure of the operators to be decided after performing the Dirac 
inversions. Rather than using pure noise in its estimators, distillation uses sources derived from the eigenvectors of 
the gauge-Laplacian, giving better coverage of relevant degrees of freedom. Increasing the number of sources in this 
scheme improves statistics faster than 1/\/N. Distillation can combined with stochastic methods, which might be 
desirable if the number of sources needed to cover the volume becomes too large. 

The distillation operator on time-slice t can be written as 

N 

a{t) = v{t)v\t) ^ a^yit) ^Y.'^L'Ht^'^Ht), (i9) 

k = l 

where the V{t) is a matrix containing the first through fc"^ eigenvectors of the lattice spatial Laplacian. The baryon 
operators involve displacements (2?i) as well as coefficients (S'aiaaas) i'^ spin space: 

XB{t) = e'^'^s^,^,^.,iv,ad)i^iv,aut^iv,aur^^{t), (20) 

where the color indices of the quark fields acted upon by the displacement operators are contracted with the antisym- 
metric tensor, and sum over spin indices. Then one can construct the two-point correlator; for example, in the case 
of proton, 

where the "baryon elemental" 

can be used for all flavors of baryon and quark masses with the same displacements on the same ensemble, and the 
"perambulator" 

rap{t',t) = VHt')M-^{t',t)V{t) (23) 

can be reused for different baryon (and meson) operators after a single inversion of the M matrix. There is a large 
factor of computational power saved by "factorizing" correlators in terms of elementals and perambulators, and the 
same elementals and perambulators can be used to contract various different correlators. 

A variations of the distillation Ref. [671 - I69] . called stochastic LapH method, has been developed to improve the 
efhciency when work on even larger volumes; a demonstration on the excited-meson and -baryon, and pion-pion 
scattering phase shift are shown in these references. 

A procedure called "pruning" can be used to reduce the number of operators when the size of a correlator matrix 
begins to get out of hand. A practical procedure is to calculate the correlators with the same source and sink 
operator (i.e., the diagonal elements of the full correlator matrix) and sort them according to some metric of their 
"individuality" . One way to systematically prune is to take the matrix of inner products of the effective masses of 
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FIG. 9: (Left) The fitted energies for tlie first-excited (top) and ground (bottom) states as functions of the number of distillation 
eigenvectors A''. (Right) A subset of G\g irrep operator correlators, grouped by their inner product. The yellow blocks indicate 
overlap > 70%. G\g is the worst case among all the irreps. 



each correlator across all time slices and sort them from there. This measures a sense in which the effective masses 
of different correlators have the same shape, implying that they also have the same excited-state content, and takes 
advantage of the black-box nature of the effective mass to eliminate ambiguity. 

The middle panel of Fig. 10 shows a subset of various Gig displacement-operator correlators sorted by the values 
of their inner products with respect to one another. Any overlap larger than 70% is marked by yellow, while the 
remaining values are depicted as magenta to blue colors. (Gig has the largest overlap amongst its operators of all the 
baryon irreps.) Notice that the inner-product matrix has a clear block structure when organized this way. Each yellow 
block consists of a set of operators that yield nearly identical effective masses as functions of time. HSC exclude many 
operators by such a selection process and pick a couple dozen operators among different quark orientations for each 
irrep. This is then refined to a smaller sub-matrix by filtering the matrix by condition numbers; a variational-method 
analysis is performed on the resulting matrix, retaining the lowest eigenstates. 

Using the "distillation" [6^ technique, HSC has been able to produce highly excited states faster than previous 
two-flavor and quenched studies. The results also has better signal-to-noise ratios than the conventional approach. 
In Ref. [55], HSC uses Nf =2-1-1 anisotropic clover lattices with pion masses of 392(4), 438(3) and 521(3) MeV at a 
fixed volume, 16^ x 128 (which gives roughly 2-fm spatial volume). 



Fig. 10 shows the results for nucleon and Delta spectroscopy with lightest pion mass around 390 MeV, classified by 



the cubic-group irreps. To improve the clarity of the graphic, we also regroup the experimental 3- and 4-star states 
according to their subduced lattice irreps. Note that multiple cubic-group states will merge into the same continuum 
state; thus, we see a few repeated states on the experimental summary plot. Out of hundreds operators, only the 
best 7 and 11 operators are chosen to calculate the full correlator matrix for the nucleon and Delta. Each principal 
correlator is fitted with a two-exponential form to get rid of possible contamination from excited states that have not 
been fully projected out by the variational method. In the end, the lowest 6 states are reported. 

In the nucleon spectrum, the first excited state from the Gig irrep is the potential Roper state and appears at 
roughly 1.9 GeV, which is heavier than its negative-parity partner, which is only around 1.7 GeV. Fig. 4 in Ref. [55] 
shows the pion mass dependence from 390 to 520 MeV, and the Roper mass decreases much faster than the S'n; it 
is likely at the physical limit, the Roper becomes lighter than its negative-parity partner but studies at lighter pion 
masses need to be performed. The higher states in Gig^u may correspond to spin-7/2 states. 

Overall, the states extracted from the lattice calculations are clustered; these states maybe degenerate in the 
continuum limit. For example, the G2 irrep corresponds to states with at least spin 5/2, which has six linearly 
independent components in the continuum limit. Each G2 state must have a partner H state with the same parity in 
order to have an interpretation as a physical state. However, since unlike irrep identity, spin is not conserved on the 
lattice, lattice discretization effects can cause the H and G2 partner states to have energy difference of O(a^). 
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The experimental Hg spectrum has a wide spread of states ranging from 1.6 to 3 GeV, while lattice data for pion 
mass around 400 MeV forms 2 clusters and ranging from 1.9 to 2 GeV; however, Fig. 4 in Ref. [55] shows that 
as pion mass decreases, the spectrum spreads out for this irrep. The Hg states correspond to five experimental 
resonances: iV(1720; 3/2+) and iV(1900; 3/2+), 7V(1680; 5/2+), iV(2000; 5/2+) and 7V(1990; 7/2+). In the Hu irrep, 
the lattice results have better determined spectrum. The four low-lying lattice states spread within 1.7-1.9 GeV may 
correspond to the experimental resonances iV(1530; 3/2"), A^(1650; 3/2") and iV(1675; 5/2"). However, the threshold 
for scattering states is near the same energy as this group of lattice states. We cannot distinguish clearly between 
single- and multi-particle contributions to the spectrum in this calculation. 

The A results from HSC at the lowest pion mass (400 MeV) is shown in the lower row of Fig. 10 As in the nucleon 
case, overall the lattice data points are higher due to the heavier pion mass used in the calculation. The extracted 
energies all decrease as the pion mass is reduced closer to the physical pion mass. In the experimental irrep figure, 
there are more degenerate states due to the cubic symmetry. 

The lowest Hg state corresponds to the ground state of A; the nearby excited state is likely to be the spin-3/2 state 
A(1600). The higher states are likely to be some combination of spin-5/2+ and 7/2+ resonances and scattering states. 
The negative-parity irrep, i?„, forms two clusters of states. The lowest two states are likely to be A(1700; 3/2~) and 
A(1930; 5/2~), while some of the higher states could correspond to spin-9/2 states, possibly mixing with the G2u 
irrep. 

The lowest two Gig states seem to correspond with the experimental resonances A(1910; 1/2+) and A(1950; 7/2+). 
However, the same spin-7/2+ could show up the spectrum of irreps Hg and G2g due to the degeneracy caused by the 
broken rotational symmetry. Similarly for the spin-9/2 state and the parity-negative irrep Giu- 

Finally, the G2 irrep starts with spin-5/2 states and then 7/2. The 3 lowest lattice states of the G2g spectrum 
might correspond to A(1905; 5/2+) and scattering states. The G2U is are less well determined, with states collapsing 
together; no strong identification with the experimental resonances can be made. 

As mentioned in the discussion of the results from the nucleon and Delta spectra, the resonance states mix with 
scattering states; thus, it is important to include multi-particle states with in the correlator matrix. Then processing 
by the variational method will help us to correctly identify the nature of the extracted excited states. Such an effort 
has been undertaken for the simplest possible scattering system, 7r7r(/ — 2). In Ref. [TOj, HSC used the distillation 
technique to study the tttt scattering process elastic phase-shift using Luscher's method. They report both the phase 
shift for isospin-2 scattering S-wave and also the D-wave channels using the same Nf = 2-1-1 anisotropic clover lattices 
with pion masses between 390 and 520 MeV. Using pairs of pions with five relative momenta (summing to zero) and 
2 smearing parameters, a 10 x 10 correlator matrix is formed and analyzed with the variational method to extract a 
spectrum of states. 

However, the above study only included the "connected" diagram, while many baryon states could mix with 
scattering states, such as and Ntttt, which include "disconnected" diagrams. In such diagrams, quarks and 
antiquarks contract within the initial state or final state. Such diagrams are notoriously difficult to calculate; not 
only do they require more sources to be calculated but also the resulting signals are much noisier. A new method for 
distillation which uses stochastic estimation of the low- lying effects of quark propagation was investigated in Ref. ^69] . 
Such an improved method resolves the problem of rapidly increasing numbers of eigenvectors as one increases the 
spatial volume, and also allows accurate determinations of disconnected diagrams. The method exploits Laplacian 
Heaviside (LapH) quark-field smearing and variance reduction is achieved using judiciously-chosen noise dilution 
projectors. They study / = 0, 1 S-wave tttt scattering using the Nf = 2 -|- 1 anisotropic clover lattices with pion mass 
as low as 240 MeV. 

The cubic-irrep operators provide many independent operators to work with on the lattice; however, we have already 
seen in the above discussion of the nucleon and Delta spectra that the multiplicity of states across different irreps can 
make them difficult to clearly identify the spin for comparison with experiments. In Ref. [56] . HSC overcomes the 
difficulty caused by the reduced symmetry of the cubic lattice by introducing a method for operator construction that 
allows for reliable identification of continuum spins; for details, see Sec. IV in Ref. Basically, we examine the 

overlap between the operator and the eigenstate in question obtained from the variational method: Z^^ = {n\Oi\0). 
The relative magnitude tells us how each operator Oi overlaps with the eigenstate n. For example. Fig. [TT] shows the 
overlaps for a set of low-lying states in the nucleon iJ„ and G2U irreps at pion mass 520 MeV and box size 2 fm. The 
overlaps for a given operator show a clear preference for overlap onto only operators of a definite spin. The assignment 
of spin must hold for states with continuum spin J subduced across multiple irreps. As can be seen, there is good 
agreement between Z values in the different irreps that are subduced from spin J after accounting for Clebsch-Gordan 
coefficients, with only small deviations. These results demonstrate that the Z values of carefully constructed subduced 
operators can be used to identify the continuum spin of states extracted for the lattices and operators shown. The 
mass values determined from fits to principal correlators in each irrep differ slightly due to the discretization effects 
and fitting variations (such as selection of the fitting intervals). Using this method, HSC determine the spectrum of 
single-particle states for spins up to 7/2, as shown in the bottom part of Fig. 11 Lighter pion masses around 400 MeV 
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FIG. 10: Nucleon (upper row, left: experiment and right: lattice with pion mass 400 MeV) and Delta (lower row, left: 
experiment and right: lattice) excited spectrum sorted according to cubic-group irrep. 



at 2 fm spatial volume are also studied in Ref . |56j , as well as a comparison of the spectrum with the non-relativistic 
quark model. 

3. Hyperon Spectra 

Using the techniques discussed in previous section, we can construct all the possible baryon interpolating operators 
that can be formed from local or quasi-local u/d and s quark fields. In Ref. [71], the masses of the lowest-lying states 
in the Gig^„ and Hg ^ representations are calculated using domain-wall valence fermions on 2-\-l fiavors of asqtad 
sea quarks (generated by the MILC collaboration[72]). The pion masses range from 300 to 700 MeV in a lattice box 
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FIG. 11: The nucleon overlap factors for a subset of operators in Hg and G2g irreps and their eigenstates (up) and extracted 
nucleon excited states using full set of operators (bottom) from Ref. . The calculation was done on the Nf = 2 + 1 anisotropic 
clover lattice at pion mass around 520 MeV. 



Flavor Gig.u 


Hg^U 


A 4 


1 


E 4 


3 


H 4 


3 


n 1 


2 



TABLE IV: The number of linearly independent interpolating operators in each irreducible representation of the cubic group 
from local-site operators 



of size 2.6 fm at lattice spacing 0.12 fm. The gauge fields are hypercubic-smeared and the source field is Gaussian- 
smeared to improve the signal. Details on the configurations can be found in Ref. [73]. The ground-state masses of 
the octet and decuplet and the SU(3) Gell-Mann-Okubo mass relation can be found in Ref. [73]. This reference has 
an extensive description of baryon-mass extrapolation (to the physical pion mass) using continuum and mixed-action 
heavy-baryon chiral perturbation theory for two and three dynamical fiavors. However, no chiral perturbation theory 
results have been published for orbitally excited hyperon resonances; in such cases, naiv e linear extrapolations in terms 
of are used. Only the local-site operators are constructed, as shown in Table the G2g,u irreps are excluded, 
since they require non-(quasi-)local operators. Figure 
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summarizes our results for the lowest-lying A and Gig^u 
(upward/downward-pointing triangles) and i?g,u (diamonds and squares) (since these have less overlap with the data 
in in Ref. [73]) and their chiral extrapolations. The leftmost points are extrapolated masses at the physical pion mass, 
and the horizontal bars are the experimental masses (if they are known). 

The lattice hyperon-mass calculations for the S, A, S and f2 are summarized on the right-hand side of Figure 12 
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FIG. 12: (left) The squared pion-mass dependence of A and fi-flavor baryons with their extrapolated values. The bar at the 
left indicates the experimental values for the corresponding spin, (right) A summary of our hyperon spectra compared to 
experiment. The lattice data from Ref. [71] are the short bars with large labels on the right. The experimental states are long 
bars with small labels on the left. 



divided into vertical columns according to their discrete lattice spin-parity irreps (Gig,u and -ffg.u), along with ex- 
perimental results obtained by subduction of continuum quantum numbers onto lattice irreps. Gi ground states 
only overlap with spin-1/2. The spin identification for H can be a bit trickier, since this irrep could match either to 
spin-3/2 or 5/2 ground states. We simply select the lowest-lying of 3/2 or 5/2 indicated in the PDG, so it could be 
either depending on which one is the ground state for a particular baryon fiavor. (We note below if the lowest H is 
not spin-3/2.) 

Although our naive extrapolation neglects contributions at ncxt-to-lcading order in chiral perturbation theory, 
several interesting patterns are seen. The better known A (with Hg being spin-5/2) and S spectra match up with 
our calculations well. Gi„ S lines up well with the S(1690), indicating that the spin-parity of this resonance could be 
l/2~; this agrees with SLAC's spin measurement. The spin assignments for the SI channel are the least known; from 
our mass pattern, we predict that 57(2250) is likely to be 3/2~ (although we cannot rule out the possibility of 5/2^), 
and 17(2380) and 17(2470) are likely spin 1/2^ and 1/2+ , respectively. 

To successfully extract even higher excited hyperon states, we would need finer lattice spacing, which would require 
massive computational resources to achieve. In Euclidean space, the excited signals exponentially decay faster than 
the ground state. One possible solution is to use an anisotropic lattice, where the temporal lattice spacing is made finer 
than the spatial ones to reduce overall costs. Such an effort have been taken by Hadron Spectroscopy Collocation and 
there are a few results reported by the HSC. Ref. [74j uses Nf = 2 anisotropic Wilson (non-O(a)-improved) lattices 
{tts/at = 3) with inverse temporal lattice spacing 5.56 GeV at a single pion mass at 416 MeV. However, only the 
lowest states of the Gig_„, G2g,u and -ffg,u of the cascade are reported. 

Higher excited states of other hyperon (S, A and S fiavors) states are followed up at low-statistics (50 to 100 
configurations) on the Nj = 2-1-1 anisotropic clover lattice in Ref. [67^. Fig. 13 shows results obtained from 2 fm box 
at pion masses around 390 MeV with eight optimal operators (from hundreds) in each symmetry channel. A detailed 
study of the Omega was reported in Ref. [3S] . Similar to the cases of the nucleon and Delta, the higher excited states 
may be subject to decay modes; thus, we need to add multiple-particle operators to the correlator matrix to separate 
them. 



IV. HEAVY FLAVOR 



Lattice QCD is now a mature field capable of providing accurate results that can be directly compared to experiment. 
In the previous sections, we showed that calculations in the light-quark sector are well established and on-going. 
However, the study of heavy quarks requires careful treatment of discretization errors. Lattice heavy quarks have 
0((toq a)") errors, where toq is the mass of the heavy quark and a is the lattice spacing. Lattice spacings for most 
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FIG. 13: Lambda (upper left), Sigma (upper right), cascade (lower left) and Omega (lower right) spectra sorted according to 
cubic-group irrep from Refs. |55l I67| . 



currently accessible dynamical ensembles are still too coarse (a~^ ~ 2 GeV) to make such systematic errors small. 
To assert better control over the discretization errors for heavy quarks on the lattice, a few specialized heavy-quark 
approaches have proven useful. For example, non-relativistic QCD (NRQCD) [75], which is an expansion of the 
lattice quark action in powers of ^J^^ , is commonly applied to bottom quarks. However, the charm-quark mass is 

not heavy enough to justify the use of NRQCD in the charm sector. Relativistic heavy-quark (RHQ) actions [7BH75] 
systematically remove 0{{mQa)"-) terms and are better suited to charm-quark calculations. Recent updates on the 
state of heavy-quark physics on the lattice can be found in several reviews [80lj87l and references therein. 

Measurements of mesons involving heavy quarks have been carried out successfully for some time and are now 
entering an era of precision. Lattice-QCD calculations of decay constants, semileptonic form factors, and neutral- 
meson mixing parameters that are important to CKM matrix elements are often quoted and considered as the 



22 



Standard Model inputs. These inputs can then be combined with experimental measurements to probe physics 
beyond the Standard Model. However, knowledge about heavy baryons is less extensive, and experimentally, there 
are still missing doubly and triply heavy baryons yet to be discovered. In this section, we emphasize recent progress 
in heavy-baryon spectroscopy from lattice QCD. 



4. Charmed Baryons 

The construction of operators in the charm sector proceeds similarly to the light-flavor ones. The interpolating 
operators used for J — 1/2 singly and doubly charmed baryons, for example, are 

■■ 73^^''' [iquCj5Qi)q's + {qfCj,Qi)qt] , (24) 

where Qu.d are the up- and down-quark fields, qs is strange-quark field and Qc is charm-quark field. 

As in the light-flavor case, people began lattice charmed-baryon calculations using the quenched approximation. 
In the early days, there were a few calculations using an 0(a)-improved light-quark action on either isotropic or 
anisotropic lattices with a single lattice spacing: Bowler et al. [88 used tree-level clover action for both light and 
heavy quarks to calculate the singly charmed baryon spectrum for spin 1/2 and 3/2. Later, Flynn et al. [89] updated 
this project with nonperturbative clover action and extended the calculation to doubly charmed baryons. Chiu 
et al. |90| used improved chiral fermion action for charm quarks and calculated both positive and negative parity 
for singly and doubly charmed baryons. Such calculations using light-quark actions to simulate heavy quarks could 
introduce large systematic errors proportional to (amg)^, which must be carefully addressed. 

One calculation has used a higher-order improved fermion action: Lewis et al. [HI] did a calculation on both doubly 
and singly charmed baryons using D234-type fermion action (which would leave a leading error of 0{a^)) for both 
light and heavy quarks but on a coarse anisotropic ensemble (with anisotropy ^ = 2). Finally, heavy-quark effective 
theory was applied to charm calculation: Mathur et al. [92] continued to use anisotropic lattices, adding two more 
lattice spacings, but changed their heavy-quark action to NRQCD, which reduces the lattice-spacing discretization 
effects. For all of these calculations, the quenched approximation remains a significant source of a systematic error 
that is difficult to estimate. 

Since new charmed baryons have been discovered in experiments, more interest has been raised to revisit charmed 
baryon calculations on the lattice. With the increasing computational resources and improved algorithms, one can use 
Nf = 2-1-1 dynamical gauge ensembles to remove the uncontrollable systematics due to the quenched approximation. 
Although more dynamical ensembles are available these days, not many charmed baryon calculations have been 
published so far. Refs. [93H95] compute the ground-state spectrum of the spin-1/2 singly and doubly charmed baryons. 
The RHQ action is used for the charm quarks and domain-wall fermions for the light valence quarks on gauge 
configurations with 2+1-flavor asqtad staggered fermions and a range of quark masses resulting in pion masses as 
light as 290 MeV at fixed lattice spacing and volume. The fermion anisotropy is nonperturbatively tuned; temporal 
and spatial clover coefficients are set separately through the tree-level tadpole improved coefficients, and two input 
bare masses for charm quarks. The results are extrapolated to the physical light-quark masses using both heavy- 
hadron chiral perturbation theory (HHChPT) as well as HHChPT-inspired polynomial extrapolations. Refs. [HHl [13 
also use 2-|-l-flavor asqtad staggered lattices but using four lattice spacings. The light-quark valence action is the 
same as the sea and the spin-3/2 charm baryons are also studied. 

All published lattice calculations on the charm baryon spectrum are summarized in Table [V] Most of the calculations 
were done in the quenched approximation with relatively heavy pion masses. A few selected mass splittings from 



dynamical charmed baryon calculations [93H97] are compared in Fig. 14 In Ref. [961 Ell > a calculation with staggered 
light quarks, an interpretation of the Fermilab action was used, defining the charm mass with the kinetic mass instead 
of the rest mass. They observe a significant lattice-spacing dependence in their calculation, which is very different 
from previous calculations [92 with quenched approximation. Their calculation covers ensembles at more than one 
lattice spacing, but their results are not as precisely determined statistically. Their uncertainties given here are merely 
the statistical ones; no systematic uncertainties are given yet. We selected a few to compare with the numbers we 



obtained in this work, as shown in Fig. 14 (For those calculations with more than one lattice spacing, we show only 
the results from the ensemble with lattice spacing closest the commonly used a « 0.12 fm among all the works.) Note 
that the quenched systematic is not included in the given errorbars. 

Finally, let us compare the doubly charmed baryons with the predictions of theoretical models, as shown in Fig. |15| 
Although the SELEX Collaboration has reported a first observation of doubly charmed baryons, searches by the 
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Nt 


Sh 


(GeV) 


L (fm) 


Bowler et al. [88] 





tree clover 1191 


2.9 


1.63 


Lewis et al. |91| 





D234 l98l 


1.8, 2.2, 2.6 


1.97 


Mathur et al. [92] 





NRQCD [19] 


1.8, 2.2 


2.64,2.1 


Flynn et al. [89] 





NP clover 


2.6 


1.82 


Chiu et al.|90l 





ODWF [99] 


2.23 


1.77 


Na et al.|96l[971 


2 + 1 


Fermilab 


2.2, 1.6, 1.3 


2.5 


Liu et al.l93ti95l 


2 + 1 


RHQ [76] 


1.6 


2.5 



TABLE V: Summary of published charmed baryon calculations from lattice QCD. Please refer to the above references and 
references within for more details. 
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FIG. 14: (Left) Comparison among charmed-baryon mass splittings of dynamical lattice calculations |93H97| . (Right) A 
summary of charmed-baryon masses (in GeV) calculated using LQCD. Two mass extractions are taken from Ref. [94]: the 
lighter (orange) points are taken from a splitting extrapolation and the darker (blue) points are taken from a direct mass 
extrapolation. 



BaBar |100| . Belle |101) and Focus |102j Collaborations have not confirmed their results. This makes it interesting 
to look back to the theory to see where the various predictions lie. We compare with a selection of other theoretical 
results, such as a recent quark- model calculation |103j . relativistic three-quark model |104], the relativistic quark 
model |105j . the heavy quark effective theory |106j and the Feynman-Hellmann theorem 107 . The lattice-computed 
mass of Sec is 3665 ± 17 ± 14 MeV, which is higher than what SELEX observed; most theoretical results suggest 
that the S^c is about 100-200 MeV higher. The flee mass prediction is 3763 ± 9 ± 44 ^35 MeV, and the overall 
theoretical expectation is for the flee to be 3650-3850 MeV. We hope that upcoming experiments will be able to 
resolve these mysteries of doubly charmed baryons. 
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FIG. 15: Comparison of theoretical predictions for doubly charmed baryons of spin 1/2. "LQCD" is the lattice QCD calculation 
taken from more precisely determined values in Refs. 194, 95 with solid error bars for the statistical error and dashed bars 
for the total error including the estimated systematic; "QM" is taken from a recent quark-model calculation [103 ; "RTQM" 
is the result of relativistic three-quark model |104| : "RQM" and "HQET" are from the relativistic quark model 105 and the 
heavy-quark effective theory [106| respectively; note that there is no error estimation done in these calculations. "FHT" is 
based on the Feynman-Hellmann theorem |107j . where rough uncertainties are estimated. 



5. Bottom Baryons 



The bottom quark is relatively heavy; thus a first-order approximation of &-quark would be to imagine a single 
static infinitely massive quark. Such an approximation introduces corrections suppressed by powers of AgcD/wf,, 
which can be included systematically in heavy-quark effective field theory. In this limit, the quark mass is removed, 
leaving only the forward hopping term in the action: 

GQ(x,t;to) = ^^ n C/4(x,t')- (25) 

The links U4 that enter this propagator can be any set of gauge covariant paths that connect neighboring points in 
time. Different choices of paths result in actions that are equivalent up to discretization errors. All of them have the 
same continuum limit. It is often accompanied by the usage of HYP-smeared links, which give significantly better 
signals for the lattice correlators. 
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In the static-quark limit, the heavy-quark propagator is reduced to a Wilson line. Thus the hadron correlator is 
quite different from the light-flavor case. For example, the static-light baryon correlation functions have the form 

Gr(x,t;^o,io) = {qfix,t)Tq']:,{x,t)ecabP''''\x,t]Xo,tQ)qf,{xo,tQ)Tq''f {xo,to)ec'a'b'), (26) 
and the static-light meson correlator is 

M{x,t;xo,to) = {q'}{x,t)P'''''{x,t;xo,to)q}'{xo,to)), (27) 

where / G {m, d, s} is the light-quark flavor index and P'^'^ is the Wilson line connecting the source and the sink which 
are separated by time t. The spin matrix F is either C75 or for scalar or vector diquarks, respectively. 

Many bottom spectroscopy lattice calculations have been performed in the meson sector. Ref. |108j calculated the 
bottomonium and bottom-light meson spectra using the relativistic heavy-quark action [78l I79j with fermion action 
parameters determined through nonperturbative tuning to eliminate systematic uncertainties; their ensembles used 
2-1-1 flavors of domain- waU fermions (DWF)[Til [TT] with hghtest pion mass 275 MeV. Ref. [109] also used DWF 
gauge ensembles but with NRQCD[75j to calculate bottom-quark quantities. A more impressive work |ll(7] done by 
the TriLat Collaboration used anisotropic 2-flavor Wilson-type lattices with renormalized anisotropy as high as 6 
and operators projected into irreducible representations of the lattice cubic group. These techniques allowed them to 
obtain very clean signals including multiple static-light excited states. 

Unlike bottomonium and B mesons, bottom baryons have not received as much attention from lattice QCD. Some 
pioneering works [88 l \92 \ 1111] were done using extrapolations with light-fermion actions or the NRQCD action in the 
quenched approximation, where the fermion loop degrees of freedom are absent. Since it is difficult to estimate the 
systematic error due to the quenched approximation, high-precision calculations cannot be achieved. Recently, more 
dynamical ensembles have become available due to the increase of the computer resources available for numerical 
research, and since the recent discoveries of the E^, S?, and ilb baryons at CDF and DO, more lattice calculations have 
emerged. Ref. [112] calculated single- and double-6 baryons with NRQCD action for the bottom quark on isotropic 
2-|-l-flavor clover ensembles, having lightest pion mass around 600 MeV. Ref. [113j also used NRQCD for bottom 
baryons but used 2-1-1 flavors of domain-wall fermions for valence/sea fermions at a single lattice spacing for spin-1/2 
and 3/2 singly to triply charmed baryons. A selection of single-5 baryons were calculated using static-quark (infinitely 
massive) action on 2-flavor chirally improved lattice Dirac operator at pion masses as light as 350 MeV in Ref. |114j . 
Ref. |115j went to a lighter pion mass (275 MeV) using 2-f 1 flavors of domain-wall fermions but used static-quark 
action to simulate the bottom quarks. Since the static-quark action tends to result in noisy signals, one needs high 
statistics to yield numbers comparable with experimental results. Ref. |116| reports static-light spin-1/2 and 3/2 
baryons on Nf = 2 twisted- Wilson fermions at a fixed lattice spacing of 0.079 fm. The lightest pion mass is 340 MeV, 
and they predict masses of negative parity states. Ref. ^117) also uses static-quark action for the bottom quarks, 
domain-wall fermion for the valence and 2-1-1 asqtad staggered fermion ensembles with lattice spacing a = 0.124 fm 
and volume about (2.5 fm)^ with lightest pion masses of 290 MeV. The number of available configurations for allowed 
this work to achieve high statistics and thus better determined resonance. An ongoing work using staggered and 
Fermilab fermion actions [7S] on 2-1-1 asqtad staggered sea and valence fermion with multiple lattice spacings was 
presented in Ref. [97] . 



A summary of results for the splittings in the bottom-hadron spectrum from Ref. |117| is shown in Table VI 
The leading-order HChPT extrapolation in terms of the mn scale is adopted as the central values, and we explicitly 
give three sources of error: statistical, extrapolation and scale. The extrapolation systematic error is estimated by 
the discrepancy between the leading-order extrapolation and the result using a form including the next-order 
term. These errors are generally about the same size as the statistical error. The error due to scale setting and 
discretization is estimated using the spread amongst the three scale-setting methods: /^r, ri and mj^. This systematic 
is quite large, up to four times the statistical error, and could be resolved by repeating the calculation on finer lattices. 
The finite-volume corrections are negligible for a 2.5-fm box, since no such correction is observed in the light-hadron 
masses (73i either. The remaining systematics, such as the effects of the AgcD/m^ corrections are omitted, since they 
are substantially smaller than the main systematics (less than one percent). 

To reduce the systematic error and because static quarks cannot produce sensible bare masses, many works present 
their results in terms of mass splittings with other nearby bottom baryons. However, not all the 2-|-l-flavor works[!)7l 
11121 11151 1117) compute the same mass splittings, so we cannot have a comprehensive comparison. Lewis et al. |112| 
used the NRQCD action to simulate their bottom quarks on CP-PACS clover lattices; Na et al.[nZ] calculated using 
the same MILC asqtad lattices as Ref. |117j but with Fermilab heavy-quark action and staggered light-quark action; 
they also included another two lattice spacings. All four calculations include the splittings S(, — A^ and Ef, — A(,; we 
calculate the additional splittings SJ^ — St, and fif, — in order to make direct comparisons. The results are shown 
in Fig. [I6] We again see good agreement amongst all lattice calculations, and mild scaling in the Na et al. results. 
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Splittmg(Stat) (Extrap) (Scale) 


Experiment 




- Bd 


71.2(2.2)(1.2)(4.2) 


87.1(0.6) 


Ab 


-Ba 


340(11)(8.1)(24) 


341.0(1.6) 




-Ba 


484.7(7.7)(6.2)(32) 


513(3) 




-Ba 


615(15)(12)(40) 


554(3) 




-Bd 


672(10)(9.1)(44) 


- 


Qb 


-Bd 


749.2(9.8) (9.0) (49) 


786(7)/886(16) 


Ai, 


-Bs 


261(10)(8.5)(19) 


253.9(1.7) 


2^ 


-Aft 


157.2(5.2)(3.4)(9.0) 


172(3) 




-At 


274(13)(13)(17) 


213(3) 




-At 


335(15)(10)(20) 






-At 


414(12)(9.5)(25) 


445(7)/545(16) 


— ij 


-Sf, 


62.6(2.6)(2.0)(3.9) 




fib 


— — () 


81.5(2.7)(3.4)(4.9) 





TABLE VI; Our consensus heavy-hadron splittings in units of MeV. The first error is statistical, the second is a systematic 
associated with the mass extrapolation and the third is a systematic associated with the setting of the lattice scale. 



Agreement with experimental values as given by the PDG is fairly good. The discrepancy in the Tih — Af, splitting 
may be due to discretization effects, as suggested by the trend visible in the Na et al. results. 

These results are in fairly good agreement with experimental results. The largest discrepancies appear for the 
Sb, although it is only about one sigma away, given the estimated systematics. We give both the DO[TlJ and new 
CDF (118j results for the fib mass. Our result (like most theory predictions) is more consistent with the lower CDF 
value. 

The fib and the yet-to-be-observed SJ^, masses are summarized along with those of other theoretical approaches in 
Fig. [171 

Finally, Meinel et al. |119j use NPQCD to simulate the bottom quark on 2 + 1 DWF and MILC asqtad ensembles. 
On each ensembles, 2 lattice spacings around 0.012 and 0.08 fm are used. The lightest sea pion ensemble is around 
300 MeV and predicts the flbbb mass to be 14.371(4)stat(ll)sys GeV. 
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FIG. 16: Comparison of mass splittings with all available 2-|-l-flavor lattice calculations of bottom baryons. The square (blue) 
points are the points extrapolated using the fi-mass reference scale; the solid error bars indicate the statistical error, and the 
dashed bars indicate the total errors (including the estimated systematic ones). The solid (red) bars indicate the experimental 
values given in the PDG, where available. For the i}t, we show both the DO result [TT] (upper-right, purple) and the CDF 
result [118] (lower-left, magenta). 
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FIG. 17: Comparison of various theoretical estimates of the Q.b (top) and Hj, (bottom) masses (and experimental values for 
the fit). "LQCD" is the lattice QCD calculation done in this work with sohd error bars for the statistical error and the dashed 
bars for the total error including the estimated systematic; "l/A^c" is taken from a combined expansion in 1/Nc, l/mq and 
SU(3) flavor symmetry breaking 120 , where the errorbar is derived solely from experimental inputs; "HQET" is the result of 
heavy-quark effective field theory (HQET) using QCD sum rules [12f| : "RQM" is from the relativistic quark model [122j : and 
"FHT" is based on the Feynman-Hellmann theorem |107| . 



SUMMARY AND OUTLOOK 



Lattice QCD has shown great potential and progress in solving the mysteries of the QCD and its resonances. Many 
experimentally well-measured quantities have been successfully reproduced. In some cases, lattice QCD can provide 
higher-precision numbers than experiments, such as the scattering length of 7r7r(/ = 2) and B-meson decay constants, 
and provide the fundamental parameters of the Standard Model, such as the strong coupling constant and quark 
masses. In nuclear physics, the H-dibaryon, a theoretical AA bound state predicted by Jaffe in the late 1970s but yet 
to be observed in experiment, has recently been reported with significant signal by two independent lattice groups. 
In addition, there are cases where experiments are impossible or extremely difficult, for example, the hyperon axial 
couplings and form factors involving initial- and final-state weak decays. Lattice QCD can provide all the necessary 
information on these quantities if provided sufficient computational resources. 

In this article, we have reviewed a few selected topics in baryon spectroscopy from lattice QCD. The mysterious 
nature of Roper state is studied through direct calculation of its mass and electromagnetic properties of its transition 
to or from the nucleon. Early studies on the Roper were mainly done in the quenched approximation, where the 
systematics are uncontrolled. In the recent years, more dynamical calculations have been reported by various groups 
(for example, HSC and CSSM) with different choices of lattice fermion actions and operator bases. The variational 
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method is used to reliably extracted the states that are nearby the Roper. Unfortunately, the latest dynamical 
calculation still have not yet observed a Roper mass smaller than its negative-parity partner. The lightest pion mass 
used so far is 156 MeV but the systematics due to finite-volume effects could be significant and need to be properly 
taken into account. At lighter pion mass and large volume, we also need to include the multiple-particle operators 
along to distinguish resonances from scattering states. Collaborations like HSC, CSSM are currently investigating 
this area, and we expects tremendous progress to be made in the near future. 

In the realm of highly radially and orbitally excited baryon states, HSC has been taking a leading role, developing 
large operators bases, special lattices (with finer temporal resolution), distillation (a novel method to speed up 
quark contractions and improve signal) and spin identification techniques (to map between cubic-group and SO (4) 
irreps). The nucleon and Delta spectra at the lightest pion mass around 390 MeV first shown in terms of cubic irreps 
classification (where one sees similar experimental results with multiple degeneracy states among irreps) and connected 
to the continuum world by examining the overlap factors with the eigcnstatc from diagonalization. As in the case 
of the Roper, we need to be able to treat all the excited sates, and being able to distinguish the non-resonance ones 
(multiple particles) is an important task. HSC has two independent studies for the mesonic case with high angular 
momentum and notoriously difficult channels involving disconnected diagrams with significant signals. Channels like 
wN and ttttN are important to baryon resonance programs and we expect more progress will soon be made, given the 
promising techniques developed by HSC. 

Hyperon excited states are also interesting, but there have been fewer calculations of them in lattice QCD. Like 
Roper, the nature of A(1405) (though not reviewed in this article) could be the first radially excited state of the A 
baryon or a five-quark KN boimd state (among other possibilities). In cascade spectroscopy, there arc resonances 
being observed experimentally whose quantum numbers remain unknown; the spin-identification carried out by the 
HSC would surely be very helpful to fill in this missing information. New experimental facilities are starting strange 
production and we expect more experimental data will be available soon to map out the unknown territory of the 
strange-baryon spectroscopy. For example, the non-zero strangeness photoproduction by CLAS at Jefferson lab is 
starting to collect high-statistics on hyperon resonances. A coherent effort by both experiment and lattice QCD would 
greatly enhance our knowledge of the hyperon spectrum. 

In the heavy-fiavor sector, the recent experimental discoveries of new charmed and bottom baryons have re-energized 
the field. They have inspired more lattice-QCD studies using variety of heavy-quark formulations, such as static ap- 
proximation, NRQCD and relativistic heavy-quark action. These new calculations have been performed on dynamical 
sea ensembles thanks to increases in the available computational resources and improved algorithms for QCD-vacuum 
generation. Many calculations are using multiple lattice-spacing ensembles to remove discretization systematic errors. 
High-precision heavy-baryon masses from lattice QCD can be achieved within the next couple years; these would 
provide useful input for experimental programs to look for the predicted resonances at specific energies. 

Much progress has been made in lattice QCD improving first-principles calculations of baryon spectroscopy in both 
light and heavy sectors. Recent developments have been discussed in this article, and we have seen the maturity of 
various techniques that will begin to resolve one of the remaining challenges: to identify the unstable resonances. 
This will involve firstly expanding the single-baryon operator basis to include multiparticle operators, sorting them 
according to their lattice symmetry and corresponding continuum spin assignment. Currently, scattering phase shifts 
for resonances and their decay modes can be computed via Luscher's method by varying the spatial volume in the 
lattice calculation. Process involving "disconnected contributions", such as 7r7r(J = 0), has shown improvement with 
the introduction of improved stochastic estimators. Process such as ttN and ttttN which are more directly relevant to 
excited baryon spectrum would be the next challenge, since the noise-to-signal on lattice QCD calculation increases 
greatly with the number of nucleons involved in the process. 

As the era of extreme-scale computing approaches, it allows lattice QCD to be able to greatly improve precision 
and attack problems that were previously considered computationally infeasible. There are multiple collaborations 
starting the gauge generation at the physical pion mass; thus, the uncertainty due to chiral extrapolation from around 
300-MeV pions to the physical mass will be nearly eliminated. This is extremely helpful to highly excited baryon 
spectroscopy, since many of the forms and low-energy constants needed for chiral extrapolation are unknown. One 
challenge is to extend the current methods for analyzing elastic processes to the treatment of inelastic decays involving 
multiple final states. This would also open avenues to investigate gluonic degrees of freedom, such as in exotic mesons 
or glueballs. There will be an intense experimental effort, such as in the GlueX experiment at the 12-GeV upgrade 
at the Jefferson Laboratory, to look for evidence of gluonic degrees of freedom in the spectrum. Calculations so 
far are limited to significantly higher pion mass than the currently available pion-mass range. In addition, due to 
the uncontrolled mixing with mesonic states, signal-to-noise is generally orders of magnitude worse than fermionic 
observables, necessitating use of the cheap quenched approximation. 

With the increasing computational resources made available to lattice QCD calculations, many obstacles have been 
overcome to reveal more of the nature QCD. Extreme-scale computing for physics has caught the attention of major 
funding agencies, and it is likely that we will see an exaflops system by 2017. For lattice QCD, this means we can 
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look forward to greater things! 
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